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1 INTRODUCTION 



Perturbative QCD (pQCD) successfully predicts inclusive energy spectra of particles in jets. To this end 
it was enough to make one step beyond the leading "Double Logarithmic Approximation" (DLA) which 
is known to overestimate soft gluon multiplication, and to describe parton cascades with account of first 
sub-leading single logarithmic (SL) effects. Essential SL corrections to DLA arise from: 

* the running coupling a s (k"j_); 

* decays of a parton into two with comparable energies, z ~ 1 (the so called "hard corrections", taken 
care of by employing exact DGLAP [ 1 1 splitting functions); 

* kinematical regions of successive parton decay angles of the same order of magnitude, 0j+i ~ 6j. The 
solution to the latter problem turned out to be extremely simple namely, the replacement of the strong 
angular ordering (AO), 8j+i <S 0«, imposed by gluon coherence in DLA , by the exact AO condition 
©j+i < ©i (see O and references therein). The corresponding approximation is known as MLLA 
(Modified Leading Logarithm Approximation) and embodies the next- to-leading correction, of order 7q, 
to the parton evolution "Hamiltonian", 70 oc ^/a~ s being the DLA multiplicity anomalous dimension 

So doing, single inclusive charged hadron spectra (dominated by pions) were found to be mathematically 
similar to that of the MLLA parton spectrum, with an overall proportionality coefficient K, ch normalizing 
partonic distributions to the ones of charged hadrons; K, ch depends neither on the jet hardness nor on the 
particle energy. This finding was interpreted as an experimental confirmation of the Local Parton-Hadron 
Duality hypothesis (LPHD) (for a review see [3 1[4| and references therein). However, in the ratio of two 
particle distribution and the product of two single particle distributions that determine the correlation, this 
non-perturbative parameter cancels. Therefore, one expects this observable to provide a more stringent 
test of parton dynamics. At the same time, it constitutes much harder a problem for the naive perturbative 
QCD (pQCD) approach. 

The correlation between two soft gluons was tackled in DLA in 0. The first realistic prediction with 
account of next-to-leading (SL) effects was derived by Fong and Webber in 1990 [6|. They obtained 
the expression for the two particle correlator in the kinematical region where both particles were close 
in energy to the maximum ("hump") of the single inclusive distribution. In [ 7 ] this pQCD result was 
compared with the OPAL e + e~ annihilation data at the Z° peak: the analytical calculations were found 
to have largely overestimated the measured correlations. 

In this paper we use the formalism of jet generating functionals [ 8 ] to derive the MLLA evolution equa- 
tions for particle correlators (two particle inclusive distributions). We then use the soft approximation 
for the energies of the two particle by neglecting terms proportional to powers of x\,x% <C 1 (x is 
the fraction of the jet energy carried away by the corresponding particle). Thus simplified, the evolu- 
tion equations can be solved iteratively and their solutions are given explicitly in terms of logarithmic 
derivatives of single particle distributions. 

This allows us to achieve two goals. First, we generalize the Fong-Webber result by extending its 
domain of application to the full kinematical range of soft particle energies. Secondly, by doing this, we 
follow the same logic as was applied in describing inclusive spectra namely, treating exactly approximate 
evolution equations. Strictly speaking, such a solution, when formally expanded, inevitably bears sub- 
sub-leading terms that exceed the accuracy with which the equations themselves were derived. This 
logic, however, was proved successful in the case of single inclusive spectra |9j, which demonstrated 
that MLLA equations, though approximate, fully take into account essential physical ingredients of 
parton cascading: energy conservation, coherence, running coupling constant. Applying the same logic 
to double inclusive distributions should help to elucidate the problem of particle correlations in QCD 
jets. 

The paper is organized as follows. 

* in section |2 we recall the formalism of jet generating functionals and their evolution equations; we 
specialize first to inclusive energy spectrum, and then to 2-particle correlations; 
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Figure 1 : Two-particle correlations and Angular Ordering 

• in section |3 we solve exactly the evolution equations in the low energy (small x) limit; how various 
corrections are estimated and controlled is specially emphasized; 

• section |4]is dedicated to correlations in a gluon jet; the equation to be solved iteratively is exhibited, 
and an estimate of the order of magnitudes of various contributions is given; 

• sectionals dedicated to correlations in a quark jet, and follows the same lines as section^ 

• in section |6] we give all numerical results, for LEP-I, Tevatron and LHC. They are commented, 
compared with Fong-Webber for OPAL, but all detailed numerical investigations concerning the size of 
various corrections is postponed, for the sake of clarity, to appendix [0 

• a conclusion summarizes this work. 

Six appendices provide all necessary theoretical demonstrations and numerical investigations. 

• in appendix and |B] we derive the exact solution of the evolution equations for the gluon and quark 
jet correlators; 

• appendix|Clis a technical complement to subsection 14. 2t 

• in appendix |D] we demonstrate the exact solution of the MLLA evolution equation for the inclusive 
spectrum and give analytic expressions for its derivatives; 

• appendix |E|is dedicated to a numerical analysis of all corrections that occur in the iterative solutions 
of the evolution equations; 

• in appendix |F| we perform a comparison between DLA and MLLA correlators. 

2 EVOLUTION EQUATIONS FOR JET GENERATING FUNCTION- 
ALS 

Consider (see Fig.Q a jet generated by a partem of type A (quark or gluon) with 4-momentum p = (j>q = 
E,p). 

A generating functional Z(E, 0; {u}) can be constructed [8] that describes the azimuth averaged parton 
content of a jet of energy E with a given opening half-angle 0; by virtue of the exact angular ordering 
(MLLA), it satisfies the following integro-differential evolution equation Q 
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d 



Z B (zp,e;{n}) Z c ((l-^)p,e;M) -Za(p, 6; {«})); (1) 



in (Q, z and (1 — z) are the energy-momentum fractions carried away by the two offspring of the 
A —> BC parton decay described by the standard one loop splitting functions 



*M(z) = C F ^, ^\z) = C F l + { \- z) , (2) 

$f(z) = T R (z 2 + (1 - z) 2 ) , = 2C A + -JL- + Z (i- Z )j , (3) 

C A = N C1 C F = (N 2 -1)/2N C , T R = l/2, (4) 

where N c is the number of colors; Z A in the integral in the r.h.s. of (Q accounts for 1-loop virtual 
corrections, which exponentiate into Sudakov form factors. 

a s (q 2 ) is the running coupling constant of QCD 

4.7T 

a s (q 2 ) = 2 — ' ^ 

QCD 

where Aqcd ~ a few hundred MeV's is the intrinsic scale of QCD, and 

/? = — f — ^ c - -n f T R ) (6) 
M 4iV c V 3 3 1 J 

is the first term in the perturbative expansion of the j3 function, nj the number of light quark flavors. 

If the radiated parton with 4-momentum k = (ko,k) is emitted with an angle G with respect to the 
direction of the jet, one has (k± is the modulus of the transverse trivector k± orthogonal to the direction 
of the jet) k± ~ \k\Q w k Q w zEQ when z -C 1 or (1 — z)EQ when z — ► 1, and a collinear cutoff 
k± > Qo is imposed. 

In (Q the symbol {u} denotes a set of probing functions u a (k) with k the 4-momentum of a secondary 
parton of type a. The jet functional is normalized to the total jet production cross section such that 

Z A (p,e-u = 1) = 1; (7) 

for vanishingly small opening angle it reduces to the probing function of the single initial parton 

Z A (p,Q^0;{u})=u A (k=p). (8) 

To obtain exclusive n-particle distributions one takes n variational derivatives of Z A over u(fcj) with ap- 
propriate particle momenta, i = 1 . . . n, and sets u = after wards; inclusive distributions are generated 
by taking variational derivatives around u = 1. 
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2.1 Inclusive particle energy spectrum 



The probability of soft gluon radiation off a color charge (moving in the z direction) has the polar angle 
dependence 

sin0d0 dsin(0/2) d9 



2(1 - cos 6) sin(0/2) ' 
therefore, we choose the angular evolution parameter to be 

r = in^^^r = ™> ; ,9) 

Qo sm(0/2) 
this choice accounts for finite angles 0(1) up to the full opening half-angle = ir, at which 

Ye=7r = In— , 
Qo 

where 2E is the center-of-mass annihilation energy of the process e + e~ — ► gg. For small angles © 
reduces to 



Tr , EQ ^ 1 d d 

Y ~ In , 6 <C 1, = , 

Qo dY dn0' 



(10) 



where is the maximal transverse momentum of a parton inside the jet with opening half-angle 0. 

To obtain the inclusive energy distribution of parton a emitted at angles smaller than with momentum 
k a , energy E a = xE in a jet A, i.e. the fragmentation function D^(x,Y), we take the variational 
derivative of Q over u a (k) and set u = 1 (which also corresponds to Z = 1) according to 



xD a A (x,Y) = E a 



5u(k a 



-Z A (k,e;{u}) 



u=l 



(11) 



where we have chosen the variables x and Y rather than k a and 0. 

Two configurations must be accounted for: B carrying away the fraction z and C the fraction (1 — z) 
of the jet energy, and the symmetric one in which the role of B and C is exchanged. Upon functional 
differentiation they give the same result, which cancels the factor 1/2. The system of coupled linear 
integro-differential equations that comes out is 



d 

dY 



xD A (x,Y) = [ dz V$ 



B A(Z)°- 

n — 



-D% ^Y + \nz)-^xD a A (x,Y) 



(12) 



We will be interested in the region of small x where fragmentation functions behave as 

xD{x) x ~ 1 p(lnx), 
with p a smooth function of In x. Introducing logarithmic parton densities 

Q = xDZ(x,Y), G = xD£(x,Y), 



(13) 



(14) 
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respectively for quark and gluon jets, we obtain from ST3 



Q{1- z)-Q) + G(z 



dy J 7T 



(G(z) - zG) + n f $£(*) (2Q(z) - d) 



(15) 



(16) 



where, for the sake of clarity, we have suppressed x and Y and only kept the dependence on the integra- 
tion variable z, e.g., 

G(z) = ?-D%( ~,Y + In z^, (17) 

such that 

G = G(1), Q = Q(1). (18) 
Some comments are in order concerning these equations. 

• We chose to express the derivative with respect to the jet opening angle 6 on the l.h.s.'s of equa- 
tions (fT3l(ll6t in terms of 

xEQ E a e 1 E 

y = Y-e = ln—— = ln-§-, £ = ln- = ln— - , (19) 

Qo Qo x E a 

instead of Y defined in ®. The variable y is convenient for imposing the collinear cutoff condition 
k± ~ xEsin 9 > Qo since, for small angles, it translates simply into y > 0; 

• to obtain d!5t one proceeds as follows. When B is a quark in (fT2l , since A is also a quark, one 
gets two contributions: the real contribution D a B=q and the virtual one — \D a A=q \ 

- in the virtual contribution, since $q(z) = &q(l — z), the sum over B cancels the factor 1/2; 

- in the real contribution, when it is a quark, it is associated with $>q(z) and, when it is a gluon, 
with &q(z); we use like above the symmetry &q(z) = &q(l — z) to only keep one of the two, 
namely <&q, at the price of changing the corresponding D(z) into D(l — z); 

• to obtain (fTol l. one goes along the following steps; now A = g and B = q or g; 

- like before, the subtraction term does not depend on B and is summed over B = q and 
B = g, with the corresponding splitting functions <& q g and Qg. In the term § 9 g , using the 
property Q 9 g (z) = &g(l — z) allows us to replace ^ J Q dz^g(z) = J z^ 9 g (z). This yields 
upon functional differentiation the — zG term in (fTBT l. For B = q,2nf flavors (nj flavors of 
quarks and nj flavors of anti-quarks) yield identical contributions, which, owing to the initial 
factor 1/2 finally yields nj; 

- concerning the real terms, & g G in dl6t comes directly from Qg^Dg in (112b . For B = q,2nf 
flavors of quarks and antiquarks contribute equally since at x <C 1 sea quarks are produced 
via gluons 3 . This is why we have multiplied Q(z) by 2nj in dl6t . 



Accompanied by a relatively small fraction 0(^/a^) of (flavor singlet) sea quark pairs, while the valence (non-singlet) 
quark distributions are suppressed as O(x). 
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Now we recall that both splitting functions $>q(z) and $ 9 g are singular at z = 0; the symmetric gluon- 
gluon splitting $>g(z) is singular at z = 1 as well. The latter singularity in ( fTBT l gets regularized by the 
factor (G(z) — zG) which vanishes at z — ► 1. This regularization can be made explicit as follows 

J dz$ 9 g (z) (G(z) - zG^j = dz& g (z) (1 — z)G(z) +z(g(z) — G 



since &g(z) = 0f (1 — z), while leaving the first term Jq dz$> g (z)(l — z)G(z) unchanged, we can rewrite 
the second 

y dz<5> 9 g {z)z(G{z) -G} = dz$|(z)(l - z)((?(l - z) - g) , 
such that, re-summing the two, (1 — z) gets factorized and one gets 



dz^(z) (g(z) -zG) = j dz^ 9 g (z)(l 



G(z)+ G(l-z)-G 



(20) 



Terms proportional to G(z) on r.h.s.'s of equations d!5l > <IT6b remain singular at z — > and produce en- 
hanced contributions due to the logarithmic integration over the region x <C z <C 1. 



Before discussing the MLLA evolution equations following from (fT3l and ill 6b . let us derive similar 
equation for two particle correlations inside one jet. 



2.2 Two parton correlations 

We study correlation between two particles with fixed energies x\ = U)\/E, X2 = oj 2 /E (x\ > x 2 ) 
emitted at arbitrary angles 0i and ©2 smaller than the jet opening angle ©. If these partons are emitted 
in a cascading process, then ©i > 02 by the AO property; see Fig.^ 



2.2.1 Equations 



Taking the second variational derivative of © with respect to u{k\) and u(k 2 ), one gets a system of 
equations for the two-particle distributions G^ and in gluon and quark jets, respectively: 



GW(z) + (qM(1-z)-QW) + G 1 (z)Q 2 (l-z)+G 2 (z)Q 1 (l-z) 



QV=jdz^& q (z) 
Gf=y dz^^g(z) (gW(z)-zGW) +G 1 (z)G 2 (1-z) 



(21) 



+ J dz^nf^z) 



2Q^\z)-G^)+2Q 1 (z)Q 2 (l-z) 



(22) 



Like before, the notations have been lightened to a maximum, such that Q( 2 ) = Q^(z = 1), G^ = 
G^(z = 1). More details about the variables on which depend are given in subsection 13.21 Now 
using dT3T > we construct the y-derivative of the product of single inclusive spectra. Symbolically, 



(QiQ2) y = Q2 I dz^ q {x) (Q 1 (l-z)-Q 1 )+G 1 (z) 



+Q l I dz^&Jx) {Q 2 (l-z)-Q 2 )+G 2 {z 
Jo 11 1 



(23) 



Subtracting this expression from d2TT i we get 
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(Q (2) - QiQ a )„ = v [g (2) (z) + (Q (2) (i - *) - Q (2) ) 

+ (G 1 (z) - Qi) (q 2 (1 - - Q2) + (G 2 (z) - Q 2 ) (Qi(l - z) - Qi 



(24) 



For the gluon jet, making use of (fT6b we analogously obtain from d22b 

(G< 2 > - G x G 2 ) y = f dz ^& g (z) \(GW(z) - zG®) + (G x {z) 

2(qW(z)-Q 1 (z)Q 2 (z) 



dz^n f & g (z) 

7T M 



G 1 )[G 2 (l-z)-G 2 
■ (G {2) - GiG 2 



+ 2Qi(«)-Gi )[2Q 2 (l-z)-G 



(25) 



The combinations on the l.h.s.'s of d2"4T) and d23t form correlation functions which vanish when particles 
1 and 2 are produced independently. They represent the combined probability of emitting particle 2 with 
•^2)2/2, •• • when particle 1 with £i,yi,... is emitted, too. This way of representing the r.h.s.'s of the 
equations is convenient for estimating the magnitude of the various terms. 



3 SOFT PARTICLE APPROXIMATION 



In the standard DGLAP region x = 0(1) (£ = O(0)), the x dependence of parton distributions is fast 
while scaling violation is small 

d e D G)Q (£,y) _ d y D G:Q {Z,y) _ 

=if) t = 0{l), - M — = ipy = 0(a s ). (26) 

U G,Q L>G,Q 

With x decreasing, the running coupling gets enhanced while the x-dependence slows down so that, in the 
kinematical region of the maximum ("hump") of the inclusive spectrum the two logarithmic derivatives 
become of the same order: 

t/jy ~ ^ = C(V57), y ~ £ ~ \Y. (27) 

This allows to significantly simplify the equations for inclusive spectra (fT5b(fT6T l and two particle corre- 
lations (1241(1251 for soft particles, xi <C 1, which determine the bulk of parton multiplicity in jets. We 
shall estimate various contributions to evolution equations in order to single out the leading and first 
sub-leading terms in J~a~ s to construct the MLLA equations. 



3.1 MLLA spectrum 

We start by recalling the logic of the MLLA analysis of the inclusive spectrum. In fact (fT5l(fT6b are iden- 
tical to the DGLAP evolution equations but for one detail: the shift In z in the variable Y characterizing 
the evolution of the jet hardness Q. Being the consequence of exact angular ordering, this modification is 
negligible, within leading log accuracy in a s Y, for energetic partons when | In z\ < \ hxx\ = 0(1). For 
soft particles, however, ignoring this effect amounts to corrections of order 0((a s In 2 x) n ) that drasti- 
cally modify the character of the parton yield in time-like jets as compared with space-like deep inelastic 
scattering (DIS) parton distributions. 

The MLLA logic consists of keeping the leading term and the first next-to-leading term in the right 
hand sides of evolution equations (IT3t<16t . Meanwhile, the combinations (Q(l — z) — Qj in (I15t and 
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(g(1 — z) — G^j in d20l produce next-to-MLLA corrections that can be omitted; indeed, in the small- x 

region the partem densities G(x) and Q(x) are smooth functions (seelT3t of lnx and we can estimate, 
say, G(l — z) — G, using (fT3l . as 

G(l - z)-G = g(^y^—,Y + ln(l- 2)) ln(l - z). 

Since ^ ~ yja~ s (see 1731. combined with a s this gives a next-to-MLLA correction C(7q ) to the r.h.s. of 
( fTBT l. Neglecting these corrections we arrive at 



Q,= r dz^& q (z)G(z), 

J x ft 

(l-z)^(z)G(z)+n f ^ g (z)(2Q(z)-G 



G y = I dz 



(28) 
(29) 



To evaluate d28t . we rewrite (see $2%) 



M(z)=C F [- + z-2 



the singularity in 1/z yields the leading (DLA) term; since G(z) is a smoothly varying function of In z 
(see (fT3l(fT4b ). the main z dependence of this non-singular part of the integrand we only slightly alter by 
replacing (z — 2)G(z) by (z — 2)G, which yields 4 



Q y = j dz ^C F (^G(z) + (z- 2)g) 



C F f 1 dz 2N c a s 



3 C F 2N c a s 



, G{z)-^-^-^^G (30) 

N c J x Z 7T W 4 A^ c 7T 



where a s = a s (In z) in the integral term while in the second, it is just a constant. To get the last term in 
d30l we used 

3 



dz(z - 2) 



(31) 



To evaluate d29l we go along similar steps. <fr q g being a regular function of z, we replace 2Q{z) — G with 
2Q - G; <& 9 g {z) also reads (see ©) 



&Jz) = 2C A 



z(l-z) 



2 + z(l-z) ; 



the singularity in 1/(1 — z) disappears, the one in 1/z we leave unchanged, and in the regular part we 
replace G(z) with G. This yields 



J dz ^ 2C A Qg(z) + (1 - z) ( - 2 + z(l - zj) G\ + n f T R (z 2 + (1 - z) 2 ) (iQ - G 



2Ca /* --<?(*) 

Z TT 



11 2 \a s 4 a s 

— Ca + ^n f T R — G + -n f T R — Q; 

6 3 / 7T 3 7T 



(32) 



the comparison of the singular leading (DLA) terms of d30t and d3*2l shows that 

„ DLA C F 



(33) 



4 since i < 1, the lower bound of integration is set to "0" in the sub-leading pieces of I28t and 1291 
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which one uses to replace Q accordingly, in the last (sub-leading) term of ( l3*2l (the corrections would 
be next-to-MLLA (seeHTl and can be neglected). This yields the MLLA equation for G where we set 

C A = N c : 



G„ 



1 dz 2N c a s 2N c a s 
— G(z) — a G 



Z 7T 



7T 



(34) 



with 



11 

12 



n f T R f 2C F \ = J_ 



11 »r 4 m ( 2C F\ 

Y N c + -n f T R (l-^ 



nt=3 



0.935. (35) 



a parametrizes "hard" corrections to soft gluon multiplication and sub-leading g — > qq splittings 5 . 

We define conveniently the integration variables z and 0' satisfying x < z < 1 and xE /Qq < &' < 6 
6 through 



£' = In — and y' = In ■ 



(36) 



The condition x < z < 1 is then equivalent to < £' < £ and xE/Qq < ©' < © is < y' < y. 
Therefore, 



^= / <1(>. 



6 d& ry 

w = L **■ 



Qo/xE 



We end up with the following system of integral equations of d3*Ul and d3*4l ) for the spectrum of one 
particle inside a quark and a gluon jet 



JO 

G(£, y) = 6(1) + £ d£' jT dy'^(£> + y') (l - a5(£' - £)) G(£', y') 



dy'll(i' + y') {G(£\ y>) - ^6(£> - £)) G(£>, y' 



(37) 



(38) 



that we write in terms of the anomalous dimension 



7o = 7o (<* s ) 



2N c a s 



7T 



(39) 



which determines the rate of multiplicity growth with energy. Indeed, using (J5J, ( fl9l and d39l one gets 

1 1 o 1 



Tol^©') 



>, , zEB'\ a z xE& 
/31n (5 ln- + — — + A 

J^QCDJ \ x Q 



W + y' + >) 



with A = ln(Qo/ Aqcd)- In particular, for z = 1 and ©' = one has 

1 1 



To 2 



P(£ + y + X) (3(Y + \y 
The DLA relation d33b can be refined to 



+ y = Y. 



(40) 



s The present formula for a differs from (47) in [ 12 1 because, there, we defined Tr = Uf/2, instead of Tr = 1/2 here. 
6 the lower bound on S' follows from the kinematical condition k± ~ xEQ' > Qo 
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,y) = §^ [l + (a - I) + «(V> 2 + ^/)) + O(7o)] G(£, y), 



(41) 



where 



G{t,y) di 
Indeed subtracting (l3lfl l and ( BTT gives 



1 dG(Ay) /2 ,. 



1 rf 2 G(l,y) 
y) t# 2 ' 



Q(£, y) - y) = §1 f o - f d V y G(l, y'); 



N r . 



(42) 



iterating twice (138b yields 



dy' 7o 2 G(£, y') = G £ + aG« + 0( 7o 2 ) = G(*, y) + a(V>| + + <% 2 ) 

which is then plugged in (|42T > to get flTT i. ip 2 + tpe? can be easily estimated from subsection 14. 21 to be 
<3(7o). In MLLA, (gTJ reduces to 



W y) = £r fi + (« - !) M*, v) + 0(7 2 )1 y). 



(43) 



3.2 MLLA correlation 



We estimate analogously the magnitude of various terms on the r.h.s. of d2"4l and d25t . Terms proportional 
to Q 2 (l—z) — Q2 and to Qi(l— z)— Q\ in the second line of d24b will produce next-to-MLLA corrections 
that we drop out. In the first line, (1 — z) — (Q^ (z) is also a smooth function of In z) will also 
produce higher order corrections that we neglect. We get 



(Q 



(2) 



QxQi)v = f dz^<S>l(z)G^(z 



(44) 



where we consider z > x% > x 2 . In the first line of (l2"5t we drop for identical reasons the term proportional 
to G 2 (l-z) - G 2 , and the term G& (z) - zG& is regularized in the same way as we did for G(z) — zG 
in (fT6b . In the second non-singular line, we use the smooth behavior of 4> q g (z) to neglect the z dependence 
in all G( 2 \ Q( 2 \ G and Q so that it factorizes and gives 



(G® - GxG 2 ) y = f dz^(l- z)M(z) G&(z) 
pi 

+ dz^ n f mz) 2(Q® - Ql Q 2 ) - ( G ® - GiG 2 ) + (2Qi - G%) (2Q 2 - G 2 ) 
Jo n L 

At the same level of approximation, we use the leading order relations 

Qi = j^ Gi, - Q,Q 2 = ^ (G( 2 ) - G 1 G 2 ) ; 



• (45) 



(46) 
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the last will be proved consistent in the following. This makes the equation for the correlation in the 
gluon jet self contained, we then get 

( G (2) _ GlG2) = f 1 dz ^1 (1 _ z) ^9( z) G (2) {z) 

J xi n 

+ [ * * nflV A (2% - 1) [(g< 2 > - (ha,) + (2% - 1) GlG2 ! 



(47) 



Like for the spectra, we isolate the singular terms 2Cf / z and 2Ca/z{\ — z) of the splitting functions <pq 
and cj)g respectively (see© and ©). We then write d44l and (147 1 as follows 



(Q (2) - QiQ 2 )„ = T ^ - 2C F [-G( 2 ) (z) + i(* - 2)G< 2 > 



(48) 



/■l r 1 

( G (2) _ GiG 2 ) y = dz^ 2C A ~G& ( z ) + (1 _ z ) C _ 2 + z(l - z)) <?( 2 ) 



/o * 



z 2 + (l-z) 2 



2 ^f - 1 



(G( 2 )-G 1 G 2 ) + (2^-l)G 1 G 2 



(49) 



which already justifies a posteriori the last equation in ( PtoT l. One then proceeds with the z integration of 
the polynomials that occur in the non-singular terms (that of d48t was already written in d3T1 i). For the 
term oc which we factorize by 2Ca, we find (see d35b for the expression of a) 



1-z) -2 + z(l-z) + 



dz 



while in the one oc G1G2 we have simply 



(50) 



Introducing 



11 n/Tfl / 2C F 



12 3JV r 



l A Nc -i nTR (l-2^ 
3 3 ; V N r 



n f =3 



allows us to express (BTT i with Ca = N c as 

a-6=^fl-^Vl-^ ^0.02, 
3iV c V N c )\ N c ) 

such that (14*81 and ( 14*91 can be easily rewritten in the form 



C F f 1 dz 2N c a s „ (2) 3C F 2N c a s ^ (2) 



v J y N c J X1 z -k 



4N C n 



(G^- Gl G 2 ) = [ ld ^ 2 ^G^\z)-a 2 -^G^ + (a-b) 2 -^ Gl G 2 

\ / y J Xl Z 7T 7T 7T 



(51) 



0.915 (52) 



(53) 



(54) 
(55) 
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Again, a s = a s (ln z) in the leading contribution while in the sub-leading ones it is a constant. We now 
introduce the following convenient variables and notations to rewrite correlation evolution equations 



, 1 , E 

In — = In — , i = 1, 2 



(56) 



y l = \n^- = \n^- = Y-£ t and r, = In ^ = £ 2 - £ 1 = Vl - y 2 > 0. (57) 
Qo Qo x 2 

The transverse momentum of parton with energy zE is k± « zEQ\. We conveniently define the in- 
tegration variables z and ©i satisfying x\ < z < 1 and ©2 < ©i < with ©2 > (©2)mm = Qo/^2 
through 



e , z x 2 EQi 
£ = in — , y = m — — 



then we write 

2/ 



r (zEQ 1 ) 



1 



2 , , X 2j E©1 

+ In h A 

x 2 



/5(ln- + ln 



7 V + y) 



(58) 



(59) 



In particular, for z = 1 and ©i = we have 

2 1 

7o " 



-1 ti 



+ Vi + V = Y. 



P(h + V2 + V + A) (3(Y + A) 

The condition xi < 2 < 1 translates into < £ < t\, while {Q 2 ) m in < ©i < © becomes < y < y 2 . 
Therefore, 



dz 



and / — — = / ay. 



a z Jo Jq /u>2 ©i Jo 

One gets finally the MLLA system of equations of d54l )( l53l for quark and gluon jets correlations 



Q®(£i,y 2 ,r l )-Q 1 (£ 1 ,y 1 )Q 2 (£ 2 ,y 2 ) = 1 f M dy^(£ + y)^--S(e-£ 1 )\G^(£,y,ri), (60) 

c J "0 



■1 ry2 



G (2) (4,y 2 ,r?) - G l (£ 1 ,y 1 )G 2 (£ 2 ,y 2 ) = d£ dy^{i + y)l-aS(£-l 1 ) G^(£,y, V ) 

Jo Jo 1 J 



U2 



+ {a-b) dyY (£ 1 +y)G(£ u y + v)G(£i+v,y)- (61) 
Jo 



In the last line of d6Tl > we have made used of d57l to write 



Gi = G(h, yi ) = G(e u y2 + v), G 2 = G(£ 2 ,y 2 ) = G{h+ V ,y 2 ). 



(62) 



The first term in d60t and d6TT) represents the DLA contribution; the terms proportional to 5 functions or 
to a, b, represent MLLA corrections, a — b appearing in doTI i and defined in ( 1531 is proportional to rif, 
positive and color suppressed. 
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4 TWO PARTICLE CORRELATION IN A GLUON JET 



4.1 Iterative solution 

Since equation doTT) for a gluon jet is self contained, it is our starting point. We define the normalized 
correlator C g by 

=C g (4, y 2 ,rj) GiG 2 , (63) 

where G\ and G2 are expressed in d62l . Substituting d63l into doTT l one gets (see appendix |A} the 
following expression for the correlator 



! _ 1 - (?i - b (-01/ + -0 2/ - jgj|]) ~ i a Xi + gg] 



1 + A + 61 + 



which is to be evaluated numerically. We have introduces the following notations and variables 



x -- 


= lnC g , 


XI 


dx 
~ d£' 


Xy 


dy' 






i>i = 


lnG l5 




1 


dGi 


Tpl,y = 


1 






dt ' 


Gi 




1p2 = 


lnG 2 , 


^2,1 


1 

G 2 


dG 2 

dl ' 


^2,y = 


1 

G 2 


dG 2 



A = 7(7 2 {i>l,ii>2,y + 1pl,yll>2,e) I 
<^2 = lo 2 (x£Xy + Xly)- 



(64) 



(65) 
(66) 
(67) 
(68) 
(69) 
(70) 



As long as C g is changing slowly with I and y, d64l can be solved iteratively. The expressions of ^ and 
ip y , as well as the numerical analysis of the other quantities are explicitly given in appendices ID.2l andlEl 
for A = (Qo = Aqcd), the so call "limiting spectrum". Consequently, (1541 will be computed in the 
same limit. 



4.2 Estimate of magnitude of various contributions 

To estimate the relative role of various terms in d64t we can make use of a simplified model for the MLLA 
spectrum in which one neglects the variation of a s , hence of 70 in d34l . It becomes, after differentiating 
with respect to I 

G ey = ^{G-aG e ). (71) 
The solution of this equation is the function for 7q = const (see appendix |C]for details) 

G(£, y) exp (270^ - «7o v) • (72) 

The subtraction term oc a in d72t accounts for hard corrections (MLLA) that shifts the position of the 
maximum of the single inclusive distribution toward larger values of I (smaller x) and partially guaran- 
tees the energy balance during soft gluons cascading (see [2][4] and references therein). The position of 
the maximum follows from d72b 

Y 

Zmax = Tri 1 + °7o)- 
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From (TTSb one gets 

= To^Ij V'j/ = ~ a ^o> V^y ~ V>« ~ fpyy = C(7o), ^ ~ V~ X = C(7o) (73) 

and the function A in d68b becomes 



A= (vS + vS)- a7o (/l + v^ 

= 2cosh( Ml -/i 2 )-a7 (e^ + e" 2 ); Mi = Ihi^. (74) 

We see that A = 0(1) and depends on the ratio of logarithmic variables i and y. One step further is 
needed before we can estimate the order of magnitude of xe, Xy an d X£y Indeed, the leading contribution 
to these quantities is obtained by taking the leading (DLA) piece of d64l . that is 

DLA , / 1 
X ~ In 1 + 



1 + A 



then, it is easy to get 

A P A 



(1 + A)(2 + A)' ^ (1 + A)(2 + A)' 
we have roughly 

Xi^Mi Xy Vy, Xiy « %; 
since ^ = /i^ = 0(7q) one gets 

X£ ~ Xy = C(7o)> ~ XiXy = O(7o), (75) 

which entails for the corrections terms <5i and <5 2 in d69b (T70b 

(5! = 0( 7o ), <5 2 = O(7o)- (76) 

The term 8\ constitutes a MLLA correction while <5 2 as well as other terms that are displayed in square 
brackets in d64l are of order 7q and are, formally speaking, beyond the MLLA accuracy. 

4.3 MLLA reduction of <|64> 

Dropping 0(7$) terms , the expression for the correlator would simplify to 

c _ 1 mlla 1 - b (Vy + rp2,i) ~ ( . 77) 



1 + A + 5 1 



4.4 C g > in the soft approximation 



C g must obviously be positive. By looking at C g > one determines the region of applicability of our 
soft approximation. Using dTTl . the condition reads 

2 + A > + (78) 

For the sake of simplicity, we employ the model d72bd73l(l74t . this gives 

2(l + cosm>i -At 2 )) > 7o(a + 6)(e' il +e^ 2 ), (79) 
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which translates into 

W- + J- >7o(a + 6). (80) 

V vi V m 

For fi,^ <^we can set y\ ~ y2 — F and, using 7q ~ 1 / /31r^ 7 , we get the condition 

v /£i'+V / ^>^-2.1, (81) 
which is satisfied as soon as > 1 (£2 > h)', so, for xi < 0.4, X2 < X\, the correlation C is positive. 

4.5 The sign of (C g - 1) 

In the region of relatively hard particles (C g — 1) becomes negative. To find out at which value of £ it 
happens, we use the simplified model and take, for simplicity, £\ = £2 = £±. 

The condition 1 = 8\ + + i>2,i)> using ( fT9l(l40l(f73l and neglecting 8\ which vanishes at £\ « £2 

reads 

\y — 77 M 4/,2 
l-& 7o . 2 J— ^ = 0^± = jL_ M g = — -4.5. (82) 

Thus in the F — > 00 limit the correlation between two equal energy partons in a gluon jet turns negative 
at a fixed value, x > x± ~ exp(4.5) = 1/90. For finite energies this energy is essentially larger; in 
particular, for Y = 5.2 (which corresponds to LEP-I energy) (I82t gives £± ~ 2.4 (x± ~ 1/11). 
For the Tevatron, let us for instance take the typical value Y = 6.0, one has £± ~ 2.6 and finally, for 
the LHC we take the typical one, Y = 7.5, one gets the corresponding £± ~ 2.8. This is confirmed 
numerically in Figs. 130 and |9] 



5 TWO PARTICLE CORRELATIONS IN A QUARK JET 

5.1 Iterative solution 

We define the normalized correlator C q by 

=C q (e 1 ,y 2 ,r,)QiQ2, (83) 

where Q\ and Q2 are expressed like in (l62l for G\ and G2. By differentiating d60l with respect to £\ and 
fj2 , one gets (see appendix IB1 



O-l 







1 - 1 (Vm + + bee] 


-[/?7o 2 ]) 


tSf gf-*l-&] 


A + 


l-f(^-[/?7o 2 ]) 


iV e Qi i " 


^-!(^-[/?7o 2 ]) 





(84) 



which is used for numerical analysis. C-JQi is computed using (|4*TI . The terms C(7q) are the one that 
can be neglected when staying at MLLA (see 15.2b . We have introduced, in addition to d65l - (f70l . the 
following notations 

(85) 
(86) 
(87) 



A = 7o \fl,W2,y + tpl,y<P2,e) , 

Si = Jq 2 Olivia + <P2,y) + 0y(<Pl,l + f2,i) 



$2 = 7o [ a l a y + °ty)i 



7 for n f =3,f3 = 0.75 
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with 



Accordingly, 

El 



ifk = lnQ k , a = lnC q . (88) 
will be computed for A = 0, the analysis of the previous functions is done in appendix 



5.2 MLLA reduction of <|S1 

Using (1351 . which entails ~ 1 - (a - § + O(7o)> reduces (TiTHi to 



C q -l 



2 + A-a(W + Vv) + [|W 
1 - a(W + - - /?7oll - *i - &] 



2 + A - + ^2,<) + [|/?7o] + + 5 2 



As demonstrated in appendix IB .21 A = A + O(7o) an(1 



C q (S 1 + h)^^C g (8 1 + S 2 ); 



and ( IH^l becomes 



f _ ~ £. 



2 + A-a(V> M + Vv) + [f/?7o]] 



(89) 



(90) 



(91) 



Would we neglect, according to d75td76b . next to MLLA terms, which amounts to dropping all 0(7g) 
corrections, d89l would simply reduce to 



mllaN c 



1 - a(ipi t i + ip 2 ,i) 



C F 2 + A - a{i) ltt + ip2,i) + Si ' 



(92) 



Furthermore, comparing {93 and ^3 and using the magnitude estimates of subsection 14.21 allows to 
make an expansion in the small 0(70) corrections 5\, ipn and tp2 i to get 



C„ — 1 MLLA N c 

Cg — 1 CF 

Nc 
C F 



1 + (b - a)(ip lt i + i/) 2 j) 



1 + A 

2 + A 



l + (6-a)(Vv + ^)(C 9 DLA )" 



(93) 



where we have consistently used the DLA expression C® LA = - — ^ . (a — 6) is given in d53b . The 

deviation of the ratio from the DLA value N c /Cf is proportional to nt, is color suppressed and numerical 
small. 
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5.3 C q > in the soft approximation 

Since we neglect NMLLA corrections and the running of a s , we can make use of d93t in order to derive 
the positivity constrain for the quark correlator. In the r.h.s. of ( 1931 we can indeed neglect the MLLA 
correction in the square brackets because it is numerically small (for instance, for 70 ~ 0.5 it is 10 -3 ). 
Therefore, C q changes sign when 

C F 5 1 
C a > 1 - -f - = - « -, 



dSOl gets therefore replaced by 



Nr. 



yi V 2/2 5 



which finally, following the same steps, gives 

5 VP 

The last inequality is satisfied as soon as l\ > 1.6 (£2 > £{)■ This condition slightly differs from that of 
the gluon correlator in 14.51 

5.4 The sign of (C q - 1) 
From <921 . C q — 1 changes sign for 



C n - 1 W 



1 -a(V>i/ + Vv) 



C F 2 + A - a(^ Xi | + ^) 



> 



(94) 



which gives the condition 



1 



This gives a formula identical to (1821 with the exchange 6 — > a; a being slightly larger than 6, we 
find now a parameter M q = 4a 2 / j3 ~ 4.66. The corresponding £± at which (C g — 1) will change 
sign is slightly higher than for gluons; for example at Y = 5.2, £± ~ 2.5 (x± ~ 1/12), Y = 6.0, 
£± ~ 2.7 (x-t ~ 1/13), Y = 7.5, £± ~ 2.9 (x-t ~ 1/16). This is confirmed numerically in figuresfJllE] 
and [10] 



6 NUMERICAL RESULTS 

In order to lighten the core of the paper, only the main lines and ideas of the calculations, and the results, 
are given here; the numerical analysis of (MLLA and NMLLA) corrections occurring in (164 1 and d84l is 
the object of appendix |EJ that we summarize in subsection 16. 3l below. We present our results as functions 

of (h + £ 2 ) and (£1 - £ 2 ). 
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6.1 The gluon jet correlator 



In order to implement the iterative solution of the first line of d64t . we define 



111 



1 + 



1 + A + [a/? 7o 2 ] 



(95) 



as the starting point of the procedure. It represent the zeroth order of the iteration for \ = lnC s . The 
terms proportional to derivatives of % m tne numerator and denominator of ( 1641 are the objects of the 
iteration and do not appear in d95l ; the parameter A depends (see d68l ) only on the logarithmic derivatives 
ipe, tp y of the inclusive spectrum G which are determined at each step, by the exact solution d!38t (I139t 
for G demonstrated in appendix |TJ] The leading piece (DLA) of (l9*5l 



nr dla , 
T g = In 



1 + 



1 



1 + A 



is the one that should be used when reducing d64l to MLLA. We have instead consistently kept sub- 
leading (MLLA and NMLLA) corrections in d95l in order to follow the same logic that proved successful 
for the single inclusive spectrum. 



6.2 The quark jet correlator 

We start now from d84l and define, like for gluons 



In < 



1 + 



C F 9 


i-i 






Cp G\ Cp G2 
N c Qi N c Q 2 


A + 




[/?7o 2 ]) 


Cf Gi , 
N c Qi "•" 


"l-!(^-[/?7o 2 ])" 


Cf G2 
N c Q 2 



(96) 



as the starting point of the iterative procedure, i.e. the zeroth order of the iteration for a = lnC q ; it again 
includes MLLA (and some NMLLA) corrections. Since the iteration concerns C q , the terms proportional 
to C g and to its derivative \t must be present in d96l l. All other functions are determined, like above, by 
the exact solution of dT38l and (f!39l for G. 

We have replaced in the denominator of d96l A with A, which amounts to neglecting O(7o) corrections, 
because the coefficient of 7(7 2 (A — A) is numerically very small; this occurs for two combined reasons: 
it is proportional to (a — 3/4) which is small, and the combination (ipi : £ y ip2,e + ip2,t,etpx,y + ip2,lyipi,e + 
ipi/iip2,y) that appears in d!17t is very small (see Fig.[T3l. Accordingly, 



^ DLA . 
T q = In 



1 + 



1 



C F 1 + A 

We can use this simplified expression for the MLLA reduction of d84l . 
6.3 The role of corrections; summary of appendix |E] 

Analysis have been done separately for a gluon and a quark jet; their conclusions are very similar. 

That tp£ and tp y , which are 0(70) should not exceed reasonable values (fixed arbitrarily to 1) provides an 
interval of reliability of our calculations; for example, at LEP-I 



2.5 < £ < 4.5 or 5 < h + £ 2 < 9, Y = 5.2. 
This interval is shifted upwards and gets larger when Y increases. 



(97) 
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Y=5.2 



Y=5.2 




0.8 1 1 1 1 1 1 1 1 1 0.9 1 1 1 1 1 1 1 1 1 

5 5.5 6 6.5 7 7.5 8 8.5 9 -2 -1.5 -1 -0.5 0.5 1 1.5 2 

h +l 2 If l 2 



Figure 2: C g for the LEP-I (Y = 7.5) inside a gluon jet as function of £\ + ^ 2 (left) and of £\ — £ 2 (right) 



T g and T q defined in d95l and d96l and their derivatives are shown to behave smoothly in the confidence 
interval Wi\ . 

The roles of all corrections 61,62, A for a gluon jet, 61, 62, A for a quark jet, have been investigated 
individually. They stay under control in ( 197 1 . While, in its center, their relative values coincide with 
what is expected from subsection 14.21 NMLLA corrections can become larger than MLLA close to 
the bounds; this could make our approximations questionable. Two cases may occur which depend on 
NMLLA corrections not included in the present frame of calculation; either they largely cancel with 
the included ones and the sum of all NMLLA corrections is (much) smaller than those of MLLA: then 
pQCD is trustable at Y = 5.2; or they do not, the confidence in our results at this energy is weak, 
despite the fast convergence of the iterative procedure which occurs thanks to the "accidental" observed 
cancellation between MLLA and those of NMLLA which are included. The steepest descent method 
[ 10 1 [ 1 1 1 , in which a better control is obtained of MLLA corrections alone, will shed some more light on 
this question. The global role of all corrections in the iterative process does not exceed 30% for Y = 5.2 
(OPAL) at the bounds of ( 197 1 ; it is generally much smaller, though never negligible. In particular, 
61 + 62 + aTg^ for gluons (or 6\ + 6 2 for quarks) sum up to 0(1O~ 2 ) at LEP energy scale (they reach 
their maximum O(10 _1 ) at the bound of the interval corresponding to the 30% evoked above). 

The role of corrections decreases when the total energy Y of the jet increases, which makes our calcula- 
tions all the more reliable. 

6.4 Results for LEP-I 

In e + e~ — > qq collisions at the Z° peak, Q = 91.2 GeV, Y = 5.2, and 70 ~ 0.5. In Fig.[2]we give the 
results for gluon jets and in Fig.|3]for quark jets. 

6.4.1 Comments 

Near the maximum of the single inclusive distribution (£1 pa i 2 ~ ¥(1 + 070)) our curves are linear 
functions of (l\ + £2) and quadratic functions of (£\ — £ 2 ), in agreement with the Fong-Webber analysis 

m 

(C q — 1) is roughly twice (C g — 1) since gluons cascade twice more than quarks (^ « 2). The difference 
is clearly observed from Fig.EJand Fig.|5](left) near the hump of the single inclusive distribution (£± + 
£2 — 7.6), that is where most of the partonic multiplication takes place. 
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5.5 6 6.5 7 7.5 8 8.5 9 -2 -1.5 -1 -0.5 0.5 1 1.5 2 

h + ! 2 h " ! 2 



Figure 3: C q for the LEP-I (Y = 7.5) inside a quark jet as function of l\ + £ 2 (left) and of i\ — £ 2 (right) 
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Y = 5.2 




- h*l2 


Y = 6.0 




-l^l 2 


Y = 7.5 



















1.55 1.6 1.65 1.7 1.75 1.8 1.85 



(l 1 + l 2 )/Y 

Figure 4: Decrease of the correlation for i\ / £ 2 at Y = 5.2, Y = 6.0 and Y" = 7.5 



In both cases, C reaches its largest value for l\ l 2 and steadily increases as a function of {l\ + i 2 ) 
(Fig.|2j left); for l\ ^ £ 2 , it increases with (£\ + l 2 ), then flattens off and decreases. 

Both C's decrease as |^i — l 2 \ becomes large (Fig. |2and|5J right). The quark's tail is steeper than the 
gluon's; for 5.9 < l\ + i 2 < 6.1, (C — 1) becomes negative when l\ — £ 2 increases; C > 1 as soon 
as £\,£ 2 > 2.75 (x\,x 2 < 0.06); this bounds is close to £ > 2.4 found in subsection !4.5l or £ > 2.5 of 
(fT50l . 

One finds the limit 

C gorq h+ ^ 2Y l. (98) 

Actually, one observes on Figs.|2j|3land|4]that a stronger statement holds. Namely, when we take the limit 
£ 2 — > Y for the softer particle, the correlator goes to 1. This is the consequence of QCD coherence. The 
softer gluon is emitted at larger angles by the total color charge of the jet and thus becomes de-correlated 
with the internal partonic structure of the jet. 
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The same phenomenon explains the flattening and the decrease of C's at i\ ^ £ 2 . 

An interesting phenomenon is the seemingly continuous increase of C g and C q at large Y for £\ pa £2 
(green curves in figs. |2]and|3]left). Like we discussed in lfT2l concerning inclusive distributions, here we 
reach a domain where a perturbative analysis cannot be trusted because of the divergence of a s . Indeed, 
when (£1 + £2) gets close to its limiting kinematical value (2Y), both y\ and 1/2 get close to 0, such 
that the corresponding a s (kf ± ) and a^kl^^) cannot but become out of control. Away from the £\ pa £2 
diagonal, taking £2 — > Y (2/2 — ► 0), we have y\ —> rj > and the emission of the harder parton still stays 
under control. 

The two limitations of our approach already pointed at in [ 12 1 are found again here: 

* x should be small enough such that our soft approximation stays valid; 

* no running coupling constant should get too large such that pQCD stays reliable. 



6.5 Comparison with the data from LEP-I 

OPAL results are given in terms of 

R(£ 1 ,£ 2 ,Y) = l + lc q (£i,£ 2 ,Y). 

In Fig.|5]we compare our prediction with the OPAL data [7 1 and the Fong-Webber curves (see subsection 
l6"6land [6|). 



6.6 Comparing with the Fong-Webber approximation 

The only pQCD analysis of two-particle correlations in jets beyond DLA was performed by Fong and 
Webber in 1990. In [6] the next-to-leading 0(70) correction, C gorq = 1 + - v /aJ + • • • , to the normalized 
two-particle correlator was calculated. This expression was derived in the region \£\ — £ 2 \/Y -C 1, 
that is when the energies of the registered particles are close to each other (and to the maximum of the 
inclusive distribution [2|[4|[ 13 1). In this approximation the correlation function is quadratic in (t\ — £2) 
and increases linearly with (£± + £2), see dlOOt . For example, if one replaces the expression of the single 
inclusive distribution distorted gaussian lTT3l (obtained in the region £ pa X- (l+a7o)) into Ml\ the MLLA 
result for a gluon jet reads 



l-(56-3^±^) 7o + 0(7o 2 



C g (lx,l 2 ,Y)r*tl + \ '- , (99) 

3 + 9^^) -(2/3 + 5a-3a^^)7o + 0(7o 2 ) 

where we have neglected the MLLA correction b\ ~ (£\ — £2) 2 \fcT s ~ near the hump of the single 
inclusive distribution {£\ pa £2 ~ ^-(1 + 070)). The Fong-Webber answer is obtained by expanding d99l 
in 70 to get [6| 



/7(FW) ^ Z 

S ~ 3 



Y 



+ 



-- (b - -a] + -3+ (b a 

3 V 3 J 9 M V 3 



Y 



70 + O(7o)- (100) 



In Fig. |6]we compare, choosing for pedagogical reasons Y = 5.2 and Y = 100, our exact solution of 
the evolution equation with the Fong-Webber predictions [6| for two particle correlations. The mismatch 
in both cases is, as seen on (llOOt . O( 7 o), and decreases for smaller values of the perturbative expansion 
parameter 70. In particular, at Y = 100, ( 7 q ~ 0.01) the exact solution doH i gets close to dlOOl ). This 
comparison is analogous in the case of a quark jet. 
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Figure 5: Correlations R between two particles produced in e + e 
and the Fong-Webber approximation 



qq compared with the OPAL data 



24 



Y=5.2 



Y=100 




5.5 6 6.5 7 7.5 8 8.5 60 80 100 120 140 160 



Figure 6: Exact C g compared with Fong-Webber's at Y = 5.2 (left) and Y = 100 (right) 



We do not perform in the present work such an expansion but keep instead the ratios d64b and ( 184b as 
exact solutions of the evolution equations. 



6.7 Predictions for Tevatron and LHC 

In hadronic high energy colliders, the nature of the jet (quark or gluon) is not determined, and one simply 
detects outgoing hadrons, which can originate from either type; one then introduces a "mixing" parameter 
to, which is to be determined experimentally, such that, the expression for two particle correlations can 
be written as a linear combination of C g and C q 

C mixed (u;l 1 ,£ 2 ,Y) = A{uj-l 1 ,t2,Y)C q {l 1 ,l 2 ,Y) + B(u;h,£ 2 ,Y)C g (£i,£ 2 ,Y), (101) 

where 



A(co;£ u £ 2 ,Y) = 



LO 


\Q(h,Y) 


Q(*2,Y)-\ 




[G(h,Y) 


G(£ 2 ,Y)\ 




r fQ(h,Y) y 
[ 1 + "{g(£ 1 ,y ) - 1 )_ 




[ (Q(£ 2 ,Y) \1 
[ 1 + "{g(£ 2 ,Y) VJ 



and 



B(iv;£ 1 ,e 2 ,Y) 



1+LO 



Q(h,Y) 
G(£i,Y) 



l+UJ 



Q(£ 2 ,Y) 
G(£ 2 ,Y) 



We plug in respectively d64b (f84b for C g and C q ; the predictions for the latter are given in Figs. f7] and [8] 
for the Tevatron, Figs.|5]andrn||for the LHC. 



6.7.1 Comments 

For both Y = 6.0 (Tevatron) and Y = 7.5 (LHC), the global behavior given in 16.4. II also holds. The 
interval corresponding to the condition C gorq > 1 is shifted toward larger values of £ (smaller x) as 
compared with the Y = 5.2 case, in agreement with the predictions of (I4.5t and <5.4b . Numerically, this 
is achieved for £ > 2.9 (£ > 3.2) at Y = 6.0 (Y = 7.5) in a gluon jet at the Tevatron (LHC). For a 
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Y=6.0 



Y=6.0 




Figure 7: C g for the Tevatron (Y = 6.0) as function of l\ + £2 (left) and of £\ — £2 (right) 
Y=6.0 Y=6.0 
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Figure 8: C q for the Tevatron (1" = 6.0) as function of l\ + £2 (left) and of £\ — £2 (right) 
Y=7.5 Y=7.5 
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Figure 9: C g for the LHC (Y = 7.5) inside a gluon jet as function of £\ + £2 (left) and of t\ — £2 (right) 
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Y=7.5 



Y=7.5 




Figure 10: C q for the LHC (Y = 7.5) inside a gluon jet as function of l x + 1 2 (left) and of l x - £ 2 (right) 



quark jet, these values become respectively £ > 3.1 (£ > 3.3) and one can check that they are close to 
the approximated ones obtained in ( 14.51 ) and < I5-4I >. 

One notices that correlations increase as the total energy (Y) increases (LHC > TeV > LEP-I). 



6.8 Asymptotic behavior of C g or q 

We display in Fig.f^the asymptotic behavior of C g and C q when Y increases. 

c ^ <»(*-l)>, 1 Y ^ <n(n-1) >, 1*, 

y < n >2 3 y < n >g 3Cf 

where n is the multiplicity inside one jet. These limits coincide with those of the DLA multiplicity 
correlator [ 14|[ 15 1. It confirms the consistency of our approach. 



7 CONCLUSION 



In this paper two particle correlations between soft partons in quark and gluon jets were considered. 

Corresponding evolution equations for parton correlators were derived in the next to leading approxima- 
tion of perturbative QCD, known as MLLA, which accounts for QCD coherence (angular ordering) on 
soft gluon multiplication, hard corrections to parton splittings and the running coupling effects. 

The MLLA equations for correlators were analyzed and solved iteratively. This allowed us to generalize 
the result previously obtained by Fong and Webber in [ 6 1 that was valid in the vicinity of the maximum 
of the single inclusive parton energy distribution ("hump"). 

In particular, we have analyzed the regions of moderately small x above which the correlation becomes 
"negative" (C — 1 < 0). This happens when suppression because of the limitation of the phase space 
takes over the positive correlation due to gluon cascading. 

Also, the correlation vanishes (C — > 1) when one of the partons becomes very soft (^ = In 1/x — > Y = 
In EQ/Qo). The reason for that is dynamical rather than kinematical: radiation of a soft gluon occurs at 
large angles which makes the radiation coherent and thus insensitive to the internal parton structure of 
the jet ensemble. 
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0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 

(ll+l 2 )/Y (l 1+ l 2 )/Y 

Figure 1 1 : Asymptotic behavior of C g and C q when Y increases 



Qualitatively, our MLLA result agrees better with available OPAL data than the Fong-Webber predic- 
tion. There remains however a significant discrepancy, markedly at very small x. In this region non- 
perturbative effects are likely to be more pronounced. They may undermine the applicability to particle 
correlations of the local parton-hadron duality considerations that were successful in translating parton 
level predictions to hadronic observations in the case of more inclusive single particle energy spectra. 

Forthcoming data from Tevatron as well as future studies at LHC should help to elucidate the problem. 
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A DERIVATION OF THE GLUON CORRELATOR C q IN & 



One differentiates — G1G2 = GiG2(C g — l) with respect to l\ and y 2 and use the evolution 
equations d38t and doTT i. 

By explicit differentiation and using the definitions (refeq:nota4bis)- (l70t one gets 



G 1 G 2 (C g -l) 



G\G<2 



+ (C g - 1) GiG 2 (^i,^2,s/ + fhtf>i,v) + G x G 2 ,ty + G 2 G 1)£ 



(102) 



the definition of x entails C g/ = xeC g , C 9i y = x y C g , C g ^ y = C g (xt y + XlXy), such that (11021 
rewrites 



G (2) - GiG 2 



C„ G\G 2 



(103) 



(Xf y + XiXy) + Xl{^l,y + i>2,y) + Xy{i>l,l + ^2,l) 

+ (C g -l) GiG 2 (ip lji ip 2 ,y + ipi t yili 2:e ) + G 1 G 2t iy + G 2 G li e y . 
By differentiating the evolution equation for the inclusive spectra d3*Sl with respect to y and i one gets 

G k ,e y = 70 (l " «(VV " /^o))^, (104) 



where one has used the definition d66td67l of to replace ^f- with Gf-ip^^, and d40l to evaluate 

1pl,l^2,y + 1pl,y1p2,t 



A 7 2 = -y3 7 *. Substituting into J103I yields 



l.h.s.m 



ty 



7 2 G 1 G 2 



(C g - 1) 2 - a ( W + 4> 2/ ) + 



+ 2a(5 1 l) +C g (5 1 + 5 2 ), (105) 



where 5\ and 5 2 are defined in (69) d70l . 

Differentiating now the r.h.s. of d6Tl with respect to y 2 and £1, one gets 



r.h.s.M 



ly 



%GiG 2 



C g (l-a (V>i/ + tp 2 ,£ - /?7o) ) ~ + (a - &) {ipi,i + - /?7o) • ( 106 ) 



Equating the expressions (11051) and (11061) for the correlation function we derive 

(C g - 1) (l + A + <5i + a ( X t + /?7o) + <5 2 ) = 1 - + ^ - /?7o) - *i - + h), (107) 

which gives (JMJ. 



B DERIVATION OF THE QUARK CORRELATOR C q IN O 

B.l Derivation of I® 



The method is the same as in appendix 1X1 one evaluates now [Q^ — Q1Q2 
First, by differentiating the evolution equation d60l . one gets 



C q - l)QiQ s 



[Q (2) - QiQ 2 ]^ = ^filc g G x G 2 (1-- (v M + ^ + xi - H)) ; (108) 



then, one explicitly differentiates 



(C q - 1)Q 1 Q 2 
C F 



and makes use of 
3 



Qk,£y = j^^G k (l - -(tpk/ ~ Pll) 



(109) 
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which comes directly from differentiating the r.h.s of d3"7t with respect to I and y; this yields 

[Q {2) -QxQi\ ly = C q Q x Q 2 



+ (C (? -l)Q 1 Q 27o 2 



-(C g - l) 7o 2 ^ [(G X Q 2 + QiG 2 ) - ^G 1 Q 2 (t/; 1/ - /? 7o 2 ) - ^QiG 2 (^ 2 , £ - /? 7o 2 ) 



(110) 



equating dl08t and dllOt gives 



C q -1 



C F U 9 


1 - § (VV + "^U + Xe~ ho) 


f ft i-c ? (i 1+ i 2 ) 


A + 




Cf G?i , 
N c Qi " + " 


i-K^-/? 7o 2 ) 


Cf G2 



(111) 



which leads (f84l . 

B.2 Expressing A, 5i and <5 2 in terms of gluon-related quantities 

All the intricacies of (lllll) lie in A, Si and J 2 defined in dSTl . which involve the quark related quantities 
(j and 99 ( l88l . In what follows, we will express them in terms of the gluon related quantities \ an d ip 



B.2.1 Expression for A 

Differentiating d4lTi with respect to I yields 



Qk,e — -Tj-Gk,. 



(112) 



then 



Qk,e f Cf „ 



Qk [N, 



1 + ( a - 7 )^*v 



+ jf°k (a - |) Vfc,« I [G k 1 - (a - |) 1 



yields 



Differentiating (|4T1 > with respect to y yields 

1 + (° ~ i)^ fc ' £ 



and, finally, 



¥>fc,j/ = tpk,y + (a - jj'fpk/y + O( 7 o)- 

Using dll4t and dll6t in A given by (f87b gives 

A W A + (a - -) (Vl/j/VV + 1p2jti>l,y + i>2Jyi>lJ + ^l,ei4>2,y) % 2 , 

which shows in particular, that A m A + C( 7 2 ). 



(113) 
(114) 

(115) 
(116) 

(117) 
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B.2.2 Expression for Si , 82 

(l87b entails Cq^Si = C q £\ipi^+(p 2 ,y) +C q> y((pi/+ip2,l)l since C q ^ and C q%y are 0(7q) and considering 
dl 16l > and dl 14b . we can approximate 

Cqllh = C q/ (^ y + lP 2 ,y) + C q>y (^ lti + W) + °(7o)> (H8) 

which needs evaluating C q e and C q>y in terms of gluonic quantities. Actually, since C q S% and C q 52 occur 
as MLLA and NMLLA corrections in (llllt . it is enough to take the leading (DLA) term of C q to estimate 
them 

differentiating then over I and y yields 

r DLA _ N c _ N c DLA 

Nr. A y Nr. 

\ + A 



Substituting dl20t . (I12U into (11181) one gets 

C g 5i =C 3 5i + O( 70 3 ). (122) 
Likewise, calculating 7 2 C g <52 needs evaluating C^ 1 /^ 4 in terms of gluonic quantities. Using dl20t one gets 

C g 5 2 =C g 5 2 + 0( 7 l). (123) 

Accordingly, C q {8\ + ^2) can be replaced by C g {5\ + S 2 ) to get the solution (111 U . This approximation 
is used to get the MLLA solution (92) of dl 1 lb - 



C DLA INSPIRED SOLUTION OF THE MLLA EVOLUTION EQUA- 
TIONS FOR THE INCLUSIVE SPECTRUM 

This appendix completes subsection 14.21 For pedagogical reasons we will estimate the solution of ( 13*81 
when neglecting the running of a s (constant^ ) (see (2][4| and references therein). We perform a 
Mellin's transformation of G(£, y) 

G(t,y) = J^^e^e-yg^v). (124) 

The contour C lies to the right of all singularities. In ( l3*8l one set the lower bounds for £ and y to 
— 00 since these integrals are vanishing when one closes the C-contour to the right. Using the Mellin's 
representation for 5(£) 



one gets 



Q{u,u) = 2— r. (126) 
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Substituting d!26t into (11241) and extracting the pole (z/ = 7 q (l/u; — a)) from the denominator of (I126t 
one gets rid of the integration over v and obtains the following representation 8 

(127) 

finally treating £ as a large variable (soft approximation x <C 1) allows us to have an estimate of dl27t 
by performing the steepest descent method; one then has 



G(£, y) U exp (2 l0 ^£y~ - a 7o 2 y) . (128) 

However, since we are interested in getting logarithmic derivatives; in this approximation we can drop 
the normalization factor of (11281) which leads to sub-leading corrections that we do not take into account 
here; we can use instead 

G(£, y) exp (2^^/Fy - a 7o 2 yj , (129) 

which is (172b . 



G(£,y)= / ^exp u£ + ^(l/u-a)y 



D EXACT SOLUTION OF THE MLLA EVOLUTION EQUATION FOR 
THE INCLUSIVE SPECTRUM 

We solve d38t by performing a Mellin's transformation of the following function ( 7 2 , and A are defined 
in (El, ©): 



F(£,y)= 1 2 (£ + y)G(£,y), 



that is, 



F(£,y) 

Substituting dl30t into d38b we obtain: 



c (2m) 



(130) 



(2m) 



dujdv 



(2mY 
f f duodv 



■ e" 1 e vy 



— a 



v UJV 

T(w,v) 



(2mJ 



where we have again replaced S(£) by its Mellin's representation (I125t . Then using the equivalence 
£ *-* -Sj, y «-> ^j, we integrate the l.h.s. by parts and obtain: 



01 



duodv 
(2m) 2 



d d , 
7r + 7T +A 



F(u,is) = P 



(2m) 2 V duJ 9u 



8 by making use of Cauchy's theorem. 
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We are finally left with the following inhomogeneous differential equation: 



\ du> dv ) v ujv v 
The variables oj and v can be changed conveniently to 

. U + V . to — V 

such that (11311) is now decoupled and can be easily solved: 



to' 2 — v' 2 ijj' — v' 



The solution of the corresponding homogeneous equation, written as a function of oj and v, is the fol- 
lowing: 

V ; pj u + s \(cj + s)v) \v + sj 

We finally obtain the exact solution of d38t given by the following Merlin's representation: 



G(£,y) = (£+y+X) —-^e^y — \ — e~ As . (132) 

J J (2m) Jo v + s\{uj + s)vJ \v + sj 

dl32l ) will be estimated using the steepest descent method in a forthcoming work that will treat two 
particles correlations at Q > Aqcd (A = Iw^Qq/Kqcd) / 0) ITOlfTTl . Substituting dl32l into doTT) 
one has the Mellin's representation inside a quark jet 

^-^im- @ - mi: M^ri^r^ 

where Jq/uju = 0(1) and the second term is the MLLA correction Jq/v = 0(70). 
D.l Limiting Spectrum, A = 

We set A = (that is Qq = Aqcd) in (I132t and change variables as follows 

u> = oj — u, s + u>t = tv/u, A = A(oj) = —— , B = a/(3 

puj 

to get (£ + y = Y is used as a variable) 



-A „ t -i 



G (£, Y) = / — x- a OY / — e QYt ( — ) t B I du u 3 ' 1 (1 + u)- A ; (133) 
the last integral of dl33t is the representation of the hypergeometric functions of the second kind (see 

M) 

rt- 1 t -B 

/ duu 3 - 1 (1 + U y A = — 2 Fi {A, B; B + 1; -t~ l ) ; 
JO 
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for ?R.B > 0, we also have 

p {n h . r . ^ _ ( Q )n 0>) n X n 

n=0 lCj « n ' 

where for example 

« = p (Q) ^ = a (a + l) ... (a + n - l) . 
Therefore (Il33t can be rewritten in the form: 

By making use of the identity fT7l : 

(l + r 1 ) 2 F! (-A + B + l,l;B + 1; -t' 1 ) = (t^) 2*1 (A, B; B + 1; -t" 1 ) , 

we split dl34t into two integrals. The solution of the second one is given by the hypergeometric function 
of the first kind CHI : 

/ — e uYt t- 1 2 Fi (-A + + 1, 1; 5 + 1; -t' 1 ) = iF (-A + 5 + 1; B + 1; -wY) . (135) 

Je2—ioo ^ % 

Taking the derivative of dl35t over (uY) we obtain: 

/ 7T- e" y * 2 Fi (-A + B + 1, 1; B + 1; -t" 1 ) =- — — ^ (-A + 5 + 1; B + 1; —uY) , 

where, 

OO / \ yi 

We finally make use of the identity [ 17 1: 

iF {-A + B + 1; B + 2; -wY) = — ^ [ iF (-A + S + 1; B + 1; -u)Y) 

" 77 ^-^ 1*1 (-A + B + 1;B + 1; -QY) 
d(—ujY) J 

to get (iF = $): 

G( "' ) =WBTT)i_, 00 + B + + 2,-010; U36) 

we can rename u) — > u; and set Y = I + y, which yields 
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G(£,y) 



+ y 



PB (B + 1) 
£ + y 



(3B (B + 1) J e _ ioo 2m 



e+ioo 

e—ioo 
e+ioo 



duJ 
2lTi 

duo 



af'S (-A + B + l,B + 2, + y)) 
e- ujy <5>(A+l,B + 2,Lu(l + y)). 



(137) 



We thus demonstrated that the integral representation (I132t is equivalent to (11361) in the limit A = 0. In 
this problem all functions are derived using (11371) . and one fixes the value of Y = ln(Q/Qo) (that is 
fixing the hardness of the process under consideration), such that each result is presented as a function 
of the energy fraction in the logarithmic scale I = ln(l/x). As demonstrated in [2 1 [ \2\, the inclusive 
spectrum can be obtained using (I136t and the result is 



G(i,Y) = 2^±K v e " BQ ^(r,y,/)J > 

1 (Y 

where the integration is performed with respect to r defined by a = - In I — — 1 ) + it, 



cosh a 



Fb(t,£,Y) 



Y a 
Z(r,£,Y) = - 



2 A 

1 — — ] sinha 



Y a 



(3 sinha 
cosh a — 



I B (2^Z(r,£,Y)), 



(3 sinha 

I B is the modified Bessel function of the first kind 



1 — — ) si nl) n 



(138) 



(139) 



D.2 Logarithmic derivatives of the spectrum, A = 

Using the expressions derived in [ 12 1 and fixing the sum I + y = Y, one gets 



d£ /J Jo 7T 



1 (1 + 2e a sinha) T B + \e a T B +i 
Y p 



and 



ay p Jo vr 



~ (1 + 2e a sinha) T B + \e a T B+ i _ 2si ^ ha ^_ 1 
Y p Y 



(140) 



(141) 



Logarithmic derivatives ip£ and ip y are then constructed according to their definition d66l)(l67l by dividing 
dl4QI > and (I14U by the inclusive spectrum (I138I) . 

Using the expression of Bessel's series, one gets 

• for £ -> 0; 



ci 



+o a . / F 
~ h ci In I — — 1 



oo, 



2a/(3+2 i-tt/2 



7r(a + 2/?) Jo 



dr (cos t) 



a//3+2 



, a \ . I a 

cos I — r — tan r sin — r 

/ V/3 



0.4097 > 0, 



-a7 + ex- 



it 



(142) 
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for £ -> Y y -> 0; 



V>f c 2 f -r - 1 ] -»■ 0, 



2 a//3+2 /-tt/ 



C-2 



7r(a + 2/?) 7 



(It ( cos r) a / 13+2 



cos ( — r ) + tan r sin ( — r 

/ V/3 



0.9218 > 0; 



~ — C2 In | — — 1 I — ► oo. 



(143) 



They are represented in Fig.^Jas functions of £ for two different values of Y (= 5.2, 15). 
D.3 Double derivatives 

In the core of this paper we also need the expression for ipj^ 
By differentiating twice dl37t with respect to £, one gets 



(144) 



Get(£,y) 



Gi(£,y) 



1 



l + y v ~ s '"" £ + y 

{£ + y)T(B) r £+io ° dtu 
+ pT(B + 3) / e _ ioo 2ri 



G(£,y) 

e-^u 2 (A 2 + 3A + 2) $ (A + 3, B + 4; + y)) . 



(145) 



Using the procedure of [ 12 1 (appendix A.2) and setting y = Y — £, the result for Gee reads 



G ee (e,Y) = y (G e (£,Y) - yG(£,Y) 
2 r ( B ) f^da e -(B-2)a 



P Jo TT 



7j2 + ^7 sinh a .Fb+i + ^2 sinh2 a 



(146) 



Likewise, for 



where 



^yy — qGw (ipy) > 



(147) 



G yy (£,y) 



-y&G(l,y) + -[G y (£,v)-G e (£,y) 



+ 



one gets 
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Figure 12: Derivatives tp£ and ifj y as functions of £ at fixed Y = 5.2 (left) and Y = 15 (right) 



r(S) /"2 da __ (B +l)a 



P Jo V 



, „ . sinh 2 a „ 1 sinh a „ 

2(5 + 1)^-^-1-^^** 



• (149) 



Finally, 

= = 70 f 1 - a fat ~ Plo)] ~ Ipilpy 

^ih ^yy an d ipe y are drawn in Fig. El of appendix IE. 1.11 as functions of £ for fixed Y. They are all 

O(7o 3 )- 

E NUMERICAL ANALYSIS OF CORRECTIONS 

In this section, we present plots for the derivatives of ijj, and tp (see d66l l (l67t and ( 1881 ). for T and 
its derivatives (see d95ld96l ). for A, Si, 62 (see (I65t - (l70b ) and the combination S c = Si + 82 + aT^, 
8 C = 81 + 82- 

E.l Gluonjet 

E.l.l ip and its derivatives 

This subsection is associated with appendices ID .21 and ID.3I . It enables in particular to visualize the 
behaviors of t/j£ and tp y when £ — > or y — > 0, as described in (11421) and d!43t . and to set the £ interval 
within which our calculation can be trusted. 

In FigElare drawn tp£ and tp y as functions of £ for two values Y = 5.2 corresponding to LEP working 
conditions, and Y = 15 corresponding to an unrealistic "high energy limit". 

ipi and (ip y ) being both 0(70), they should not exceed a "reasonable value"; setting this value to 1, 
< 1 and \tp y \ < 1 set, for Y = 5.2, a confidence interval 

2.5 < £ < 4.5. (150) 
In the high energy limit Y = 15, this interval becomes, 4.5 < £ < 13, in agreement with 14.51 
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Figure 13: Double derivatives ^„ and as functions of £ at fixed Y 



E.1.2 A(£i,£ 2 ,F) 

A has been defined in d68b . in which ipx^ and ipi y are functions of £\ and Y, ip2,e and ^2,2/ are functions 
of £ 2 and V. 

Studying the limits £ — > and £ — ► Y" of subsection lD.2l 



for 



1, «2 



such that 



Y" one gets (using the results oflDT 



A ~ -4 



Y-h, Y-i 2 Y 
— 7, hi — 1 



*2, Y 
— In — 



A^^T°0, A 



+00. 



(151) 



(152) 



for £1, £2 — > one gets (according to|D 
A 



-070 



0/1 1 \ A y- 



-00. 



(153) 



In Fig. [H (left) A is plotted as a function of £ x + £2 for three different values of l x - £2 (0.1, 0.5, 1.0); 
the condition (11501) translates into 

5.0 < £1 +£ 2 < 9.0; (154) 

on Fig.EKrighf) the asymptotic limit A — ► 2 for very large Y clearly appears (we have taken £\ — £2 = 
0.1); it is actually its DLA value \2\; this is not surprising since, in the high energy limit 70 becomes 
very small and sub-leading corrections (hard corrections and running coupling effects) get suppressed. 



E.1.3 T g and its derivatives 



Fig. [T31 exhibits the smooth behavior of exp(T 9 ) as a function of {£\ + £2) in the whole range of 
applicability of our approximation (we have chosen the same values of (£\ — £2) as for Fig. [TUl. and 
as a function of {£\ — £2) for three values of [£\ + £2) (6.0, 7.0, 8.0). So, the iterative procedure is safe 
and corrections stay under control. 
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Figure 15: exp (T g ) as a function of l\ + £2 (left) and, l\ — £2 (right) for Y = 5.2 



39 



Y=5.2 





T gJ ,0<l 1 


■l 2 <0.1 




T < I-, 
g,y u % '1 


I - 1 




T fti.y 0<l i 


-I 2 <0.1 











5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 



Y=5.2 



0.3 


T g 0.9 < I, 


- 1 2 < 1.1 






\ \ T 9,r°- 9<I 1 


- 1 2 < 1.1 




0.2 


- \ \ T.| V ,0.9<I 1 


- 1 2 < 1.1 




0.1 

















0.1 








0.2 









5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 



l 1+ l 2 



Figure 16: T ' g £, T 9 y and T g ,£y as functions of I1 + I2 for Y = 5.2, £± — £2 = 0.1 (left) and £\ — £2 = 1-0 
(right) 



Fig- El displays the derivatives of T g . dl55t . (I156t and (I157t have been plotted at Y = 5.2, for (£\ — 
£2) = 0.1 (left) and {£\ — £2) = 1-0 (right). The size and shape of these corrections agree with our 
expectations (T g/ = T 9iV = 0(7^), T g/y = 0(j$)). 

For explicit calculations, we have used 

T [l-b (V^+Vy- hi)] (A^-q/3 2 7 4 ) ^(^i,»+V'2,«+/? 2 7o 4 ) 

9/ (l + A+a/3 7 2 ) [2+A-6(^i / +Vv-/37 2 )] 2+A-6 ' 

T [1-6 (V^^y-/^ 2 )] (A,-a/? 2 7 4 ) b (Vy.+Vy ^/^p 4 ) 

"• y (l + A+a/3 7 2 ) [2+A-6(^i / +Vv-/37o 2 )] 2+A-6 ' 

T 9,ey=—gf-^ ( 157 ) 



where 



A* = 7o 2 bPlje^y + i>l,£^2,y£ + ^2,££^l,y + ^2,£^l,y i] + ^7o A, 

Aj/ = 7o~ 2 bPhly^y + lpl,£lp2,yy + ^2,£y^l,y + 1p2,£tpl,yy] + hi A- (158) 



For the expressions of = an d ^yn, the reader is directed to lD.3l ( 1157b has been computed 

numerically (its analytical expression is too heavy to be easily manipulated). 



E.1.4 8 U 8 2 ,8 C 

5\ and 82 are defined in (l65ld70l . We also define 

5c = 81 + 5 2 + aTi, (159) 
which appears in the numerator of the first line of (i64t, 

Fig.[l7]displays the behavior of 81, 8 2 and 81 + 8 2 at F = 5.2 for 4 - £ 2 = 0.1 and ^ - £ 2 = 1.0. We 
recall that these curves can only be reasonably trusted in the interval (U54t . 

Though \8\\ = O(7o) (MLLA) should be numerically larger than \8 2 \ = C(7q) (NMLLA), it turns out 
that for relatively large 70 ~ 0.5 (Y=5.2), |<5i| ~ \8 2 \, and that strong cancellations occur in their sum. 
As 70 decreases (or Y increases) 1 5\ | 3> 1 62 | , in agreement with the perturbative expansion conditions. 
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Figure 17: b x , 5 2 and ft + J 2 as functions of 4 + ^ 2 for £ 1 -£ 2 = Q.l (left) and 4 - £ 2 = 1.0 (right) 



0<l 1 -l 2 <0.1 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
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Figure 18: S c as a function of {l x + £ 2 )/Y for ti - £2 = 0.1 (left) and £ x - £ 2 = 1.0 (right) 



In Fig. [23 we represent 5 C for different values of Y; it shows how the sum of corrections (MLLA and 
NMLLA) stay under control in the confidence interval dl54t . For Y = 5.2 one reaches a regime where it 
becomes slightly larger than 0.1 away from the region x\ « x 2 (see upper curve on the right of Fig.[T8b 
but still, since 1 (which is the leading term in the numerator of doH l) S> 0.1, our approximation can be 
trusted. 



E.1.5 The global role of corrections in the iterative procedure 

Fig-EDshows the role of 5 C on the correlation function: we represent the bare function exp T g (seel95l 
as in Fig.[Bl together with {64}. For (£i - £ 2 ) = 0.1 {£1 m £ 2 ) and (£1 - £ 2 ) = 1.0, it is shown how 5 C 
modifies the shape and size of exp T g . When l\ ^ £ 2 {{£\ — £ 2 ) = 1.0), 5 C decreases the correlations. 
They are also represented as a function of (£\ — £ 2 ) when {£\ + £ 2 ) is fixed ( to 6.0 and 7.0). The increase 
of S c as one goes away from the diagonal £\^£ 2 (see Fig.^Jfor {£\ — £ 2 ) = 1.0) explain the difference 
between the green and blue curves; this substantially modifies the tail of the correlations. 

When Y gets larger, the role of 5 C decreases: at Y = 7.5 (LHC conditions) the difference between the 
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Figure 19: C g (blue) compared with exp T g (green) 



two curves becomes negligible. 

E.2 Quark jet 

E.2.1 ip and its derivatives 

Fig. |20] displays the derivatives ipe and tp y together with those tp£ and tp y for the gluon jet, at Y = 5.2. 
There sizes and shapes are the same since the logarithmic derivatives of the single inclusive distributions 
inside a gluon or a quark jet only depend on their shapes (the normalizations cancel in the ratio), which 
is the same in both cases. The mismatch at small t between ipp and ipg stems from the behavior of ip££ 

Tpi£ —oo. Therefore, in the interval of applicability of the soft approximation dll4t and dll6t can 
be approximated by tp£ and ip y respectively. 

E.2.2 A(£ 1 ,£ 2 ,Y) 

The last statement in lE.2.1l numericallv supports the approximation (II 17t . that is 

A« A + O( 7o 2 ). 
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Figure 20: Derivatives ip g i and ip gy as functions of £ at fixed Y = 5.2 (left), compared with ipg and ip, 

Y=5.2 Y=5.2 




4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 -2 -1.5 -1 -0.5 0.5 1 1.5 2 



h + ! 2 h - ! 2 

Figure 21: exp (T g ) as a function of £\ + ^2 (left) and, £\ — £2 (right) for Y = 5.2 



We get rid of the heavy 0(7q) factor in dll7t to ease our numerical calculations. Hence, the behavior of 
A is already given in Fig. EI 

E.2.3 T q and its derivatives 

The smooth behavior of expT g is displayed in Fig. |^ as a function of the sum (l\ + £2) for fixed 
(£1 — I2) and vice versa. The normalization of (exp T q — 1) is roughly twice larger (x ^ 2) than 
that of (exp T g — 1). We then consider derivatives of this expression to get the corresponding iterative 
corrections shown in Fig. [22] The behavior of Tg^O^)), Yq i3/ (0(7g)) and T g,£j/(C(7o)) i s i n good 
agreement with our expectations as far as the order of magnitude and the normalization are concerned 
(see also Fig.[l6t 9 . 

9 it is also important to remark that T q ,i, T q ,e, Yq.ey are x ^-^g,e, ^g,e, Tg,ty ■ 
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Y=5.2 



Y=5.2 





Figure 22: Corrections T q> £, T qy and „ as functions of l\ + l 2 f° r ^1 — ^2 = 0-1 (left) and l\ — t 2 = 
1.0 (right) at Y = 5.2 



E.2.4 <Ji, <5 2 and <5 C 

We define 

lc = Si + 5 2 

as it appears in both the numerator and denominator of d84t . In Fig.|23]are displayed 5\, 5 2 and their sum 
8 C as functions of the sum {h + £ 2 ) at fixed [l x - t 2 ) (h -£2 = 0.1, left) {t x -£ 2 = 1.0, right). 

At Y = 5.2, which corresponds to 70 ~ 0.5, the relative magnitude of 8% and 5 2 is inverted 10 with 
respect to what is expected from respectively MLLA and NMLLA corrections (see subsection I4.2t . This 
is the only hint that, at this energy, the expansion should be pushed to include all NMLLA corrections to 
be reliable. 

Large cancellations are, like for gluons, seen to occur in S c , making the sum of corrections quite small. 
In order to study the behavior of 5 C as Y increases, it is enough to look at Fig. |24] where we compare 8 C 
at Y = 5.2, 6.0, 7.5. 

E.2.5 Global role of corrections in the iterative procedure 

It is displayed in Fig.|25] 5 C does not affect exp T q near the main diagonal (t\ = £ 2 ), but it does far from 
it. We find the same behavior as in the case of a gluon jet. 



F COMPARING DLA AND MLLA CORRELATIONS 

In Fig. |26]we compare the quark correlator at DLA and MLLA. The large gap between the two curves 
accounts for the energy balance that is partially restored in MLLA by introducing hard corrections in 
the partonic evolution equations (terms oc a, b and |); the DLA curve is obtained by setting a, b and | 
to zero in ( E)4l and (l%4l ; C q is a practically constant function of t\ + £ 2 in almost the whole range, and 
decreases when l\ + l 2 — > 2Y by the running of a s . The MLLA increase of C q with l\ + t 2 follows 
from energy conservation. Similar results are obtained for C g . 



10 it has been numerically investigated that the expected relative order of magnitude of 5i and 82 is recovered for Y > 8.0 
(this value can be eventually reached at LHC). 
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Figure 23: Ji, <5 2 and <5i + <5 2 as functions of l x + £ 2 for i\ - £ 2 = 0.1 (left) and 4 - 4 = 1.0 (right) at 
Y = 5.2 




Figure 24: 5 C as a function of {l x + ^ 2 )/Y for 4 - £ 2 = 0.1 (left) and l x - £ 2 = 1.0 (right) at 
Y = 5.2, 6.0, 7.5 
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Figure 26: Comparing DLA and MLLA correlations 



46 



List of Figures 



1 


Two-particle correlations and Angular Ordering! 


. . . . 4 


2 


Cn for the LEP-I (Y = 7.5) inside a eluon iet as function of £1 + £o (left) and of L 


- £0 (riehd 21 


3 


Cn for the LEP-I (Y = 7.5) inside a quark iet as function of £i + £9 (left) and of £1 


- £0 (riehd 22 


4 


Decrease of the correlation for £a £0 at Y = 5.2. V = 6.0 and Y = 7.3 .... 


22 




5 


Correlations R between two particles produced in e + e — ► qq compared with the OPAL data and the Fong-Webb 


6 


Exact C n compared with Fong-Webber's at Y = 5.2 (left) and Y = 100 (right)! . . 


25 


7 


C„ for the Tevatron (Y = 6.0) as function of L 4- £0 Heft) and of L - £0 (right)! . 


26 


8 


C„ for the Tevatron (Y = 6.0) as function of L 4- £0 (left) and of £a - £0 (rieht)! . 


. . . . 26 


9 


C„ for the LHC (Y = 7.5) inside a gluon jet as function of £a + £0 (left) and of £a - 


- £0 (right)l 26 


10 


C for the LHC (Y = 7.5) inside a eluon iet as function of £1 + £0 (left) and of £1 - 


- £0 (riehd 27 


11 


Asymptotic behavior of C„ and C q when Y increases] 


28 


12 


Derivatives ibe and if) v as functions of £ at fixed Y = 5.2 (left) and Y = 15 (riehd 


. . . . 37 



13 Double derivatives tbt 0. ibt> v and ib vv as functions of £ at fixed Yl 38 



14 A as a function of li + £0 for Y = 5.2 (left) and its asymptotic behavior (rieht. It - £0 = 0.1)1 39 



1 5 exn (T J as a function of & 4- lo (left) and. & - lo (right) for Y = 5.2l 39 



1 6 xZZ Ty and T = , v as functions of £, 4- lo for Y = 5.2. L - £0 = 0.1 (left) and £ - £0 = 1 .0 (right! 40 



17 and 4- ^9 as functions of & 4- & for £, - £0 = OA (left) and li - lo = 1 .0 (right)! 41 



1 8 & as a function of (li 4- lo)/Y for L - £•, = 0.1 (left) and It - l» = 1 .0 (right)! .... 41 



19 C„ (blue) compared with exp T„ (green)! 42 



20 Derivatives u>„ v and <j>„ v as functions of £ at fixed Y = 5.2 (left), compared with ibi and q/u 43 



21 exp (T„) as a function of li 4- £0 (left) and, li - £0 (rieht) for Y = 5.j 43 



22 Corrections T„ 0. T„ and T„ e v as functions of li + £9 for £ i — I9 = 0.1 (left) and li — £0 = 1.0 (rieht) at Y 



23 S • S and ft 4- £9 as functions of jj + lo for li - lo = 0.1 (left) and & - l? = 1.0 (rieht) at Y = 5~1 45 



24 S as a function of fli + £0) IY for It - £0 = 0.1 (left) and £i - £0 = 1.0 (rieht) at Y = 5.2. 6.0. zjj 45 



25 C„ (blue) compared with exp T 7 (ereen)! 46 



26 Comparing DLA and MLLA correlations! 46 



47 



References 

[1] V.N. Gribov and L.N. Lipatov: Sov. J. Nucl. Phys. 15 (1972) 438 and 675; 
L.N. Lipatov: Sov. J. Nucl. Phys. 20 (1975) 94; 

A.P Bukhvostov, L.N. Lipatov and N.P Popov: Sov. J. Nucl. Phys. 20 (1975) 286; 
G. Altarelli and G. Parisi: Nucl. Phys. B 126 (1977) 298; 
Yu.L. Dokshitzer: Sov. Phys. JETP 46 (1977) 641. 

[2] Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller and S.I. Troyan: Basics of Perturbative QCD, Editions 
Frontieres, Paris, 1991. 

[3] Yu.L. Dokshitzer, V.A. Khoze, S.I. Troyan and A.H. Mueller: Rev. Mod. Phys. 60 (1988) 373-388. 

[4] V.A. Khoze and W. Ochs: Int. J. Mod. Phys. A12 (1997) 2949. 

[5] Yu.L. Dokshitzer, VS. Fadin and V.A. Khoze: Z. Phys C18 (1983) 37. 

[6] CP Fong and B.R. Webber: Phys. Lett. B 241 (1990) 255. 

[7] OPAL Collab.: Phys. Lett. B 287 (1992) 401. 

[8] K. Konishi, A. Ukawa, and G. Veneziano: Nucl. Phys. B 157 (1979) 45. 

[9] OPAL Collab., M.Z. Akrawy et al.: Phys. Lett. B 247 (1990) 617; 
TASSO Collab., W. Braunschweig et al.: Z. Phys. C 47 (1990) 198; 
L3 Collab., B. Adeva et al.: L3 Preprint 025 (1991). 

[10] R. Perez-Ramos: "Two particle correlations inside one jet at "Modified Leading logarithmic Ap- 
proximation" of Quantum Chromodynamics; II: steepest descent evaluation", in preparation. 

[11] R. Perez-Ramos: These de Doctorat, Universite Denis Diderot - Paris 7 (2006). 

[12] R. Perez-Ramos & B. Machet: "MLLA inclusive hadronic distributions inside one jet at high energy 
colliders", hep-ph/0 512236| JHEP 04 (2006) 043. 

[13] CP Fong and B.R. Webber: Phys. Lett. B 229 (1989) 289. 

[14] E.D. Malaza and B.R. Webber: Nucl. Phys. B 267 (1986) 702. 

[15] Yu.L. Dokshitzer, VS. Fadin and V.A. Khoze: Phys. Lett. B 115 (1982) 242. 

[16] I.S. Gradshteyn and I.M. Ryzhik: Table of Integrals, Series, and Products, Academic Press (New 
York and London) 1965. 

[17] L.J. Slater, D. Lit., Ph.D. Thesis: Confluent Hypergeometric Functions, Cambridge at the Univer- 
sity Press, 1960. 



48 



